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Through the modelling of the Spectral Energy Distribution of blazars we can infer the physical parameters required 
to originate the flux we see. Then we can estimate the power of blazar jets in the form of matter and fields. These 
estimate are rather robust for all classes of blazars, although they are in part dependent of the chosen model (i.e. 
leptonic rather than adronic). The indication is that, in almost all cases, the carried Poynting flux is not dominant, 
while protons should carry most of the power. In emission line blazars the jet has a comparable, and often larger, 
power that the luminosity of the accretion disk. This is even more true for line-less BL Lacs. If the jet is structured at 
the sub-pc scale, with a fast spine surrounded by a slower layer, then one component sees the radiation of the other 
boosted, and this interplay enhances the Inverse Compton flux of both. Since the layer emission is less beamed, it can 
be seen also at large viewing angles, making radio-galaxies very interesting GLAST candidates. Such structures need 
not be stable components, and can form and disappear rapidly. Ultrafast TeV variability is challenging all existing 
models, suggesting that at least parts of the jets are moving with large bulk Lorentz factors and at extremely small 
viewing angles. However, these fast "bullets" are not necessarily challenging our main ideas about the energetics and 
the composition of the bulk of the jet. 
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1. Introduction 

EGRET, onboard CGRO, and the ground based 
Cherenkov telescopes showed that blazars are the 
most powerful high energy extragalactic emitters, 
and allowed to know the bolometric luminosity of 
these objects. Now, just after the launch of AGILE 
and just before GLAST, we are preparing for (and 
already tasting) the possibility to have simultaneous 
data in the optical, X-rays and the GeV bands (and 
possibly the TeV one). This will be possible mainly 
by Swift: its rapid slew and flexible scheduling will 
ensure good quality optical and X-ray data while the 
7-ray observations are still ongoing. What was an ex- 
ception in the EGRET era, will be routine. This is 
therefore a risky and at the same time healthy time 
to put forward new ideas concerning blazars, that 
GLAST can falsify. In this contribution, I will try 
to summarize how the power of jets can be derived, 
and what inferences can we draw from that. I will dis- 
cuss the ultra-fast TeV variability recently observed 
in PKS 2155-304, arguing that, contrary to previ- 
ous claims, it is unlikely that the jet of this source 
is magnetically dominated. Furthermore, I will point 
out that the fact that the TeV emitter BL Lac ob- 
jects are also the least powerful blazars opens up the 
possibility to slow down their jets by the Compton 
rocket effect. 



2. From the SED to the jet power 

Fig. 1 shows three states of 3C 454.3, fitted with a 
simple homogeneous, one zone, synchrotron inverse 
Compton model [18]. The main difference among the 
three states is the bulk Lorentz factor F, assumed 
to be a function of where the jet produces most 
of the flux. This model, proposed and discussed by 
Katarzynski & Ghisellini [23] , assumes that when the 
jet dissipates closer in, it has a smaller F than when 
it dissipates at larger distances from the black hole. 
This has a great impact on the ratio between the 
inverse Compton to synchrotron peak levels Lc/Ls- 
for small F and more compact dissipation regions, 
it is likely that the magnetic held is larger, and the 
external photons, thought to originate in the broad 
line region (BLR) are seen less boosted, resulting in 
a relatively small L^jL^. Viceversa for larger F. As 
can be seen, the resulting spectral energy distribu- 
tion (SED) can change dramatically in selected fre- 
quency ranges, without the need to a large change of 
the overall jet power. This model can be contrasted 
by the recent idea by Sikora et al. [29], in which the 
2007 high 7-ray state is attributed to a dissipation 
region even beyond the BLR. At these distances, the 
magnetic field is small, and the external radiation is 
produced by a ~10 pc scale dusty torus. Variabil- 
ity of the 7-ray flux is clearly a diagnostic between 
the two possibilities. Fig. 1 shows also that when we 
see a very hard X-ray spectrum we have the best 
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Fig. 1. SEDs of 3C 454.3 at different epochs. Top panel: 
the SED in 2000, corresponding to the BeppoSAX observa- 
tions discussed in [33]. The other points are not simultaneous 
(see [26] and references therein). Mid panel; the SED dur- 
ing the huge optical flare in 2005, as described in [20, 26]. 
Bottom panel: the SED on July 26, 2007, as observed by AG- 
ILE [6] , and Swift. The optical flux in the R band comes from 
the Tuorla observatory. The optical and X— ray fluxes are cor- 
rected for Galactic extinction [Ay = 0.355). The solid and 
dashed lines correspond to our modelling. The dotted line is 
the contribution from the accretion disk (assumed to be a sim- 
ple black-body). From [18]. 

chance to measure the low energy cut-ofF of the elec- 
tron distribution, emitting, by the external Compton 
mechanism, in the X-ray range. In particular, for 3C 
454.3, the electron distribution is required to extend 
down to 7min ~ 1- This is the energy region where 
most of the electrons are, and knowing 7min is crucial 
to determine the power that the jets carries in the 
form of bulk motion of particles. 

We do not have the same information for BL Lac 
objects, that can be fitted by the SSC model without 
the need of an external component (but see below 
from the alternative spine-layer scenario). On the 
other hand, for these sources it is less important to 
know 7min, since the average energy of their electrons 
is higher, approaching the proton rest mass energy. 



The kinetic power of the jet in this case depends less 
on the total number of electrons (and protons) and 
more on the electron mean energy. 
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Fig. 2. The power of blazar jets in bulk motion of protons 
(Lp, assuming one proton per emitting electrons), relativistic 
electrons (Le), magnetic field (Lb) a^nd ™ radiation L^. The 
star and the triangle correspond to the values of 2155—304 
during the TeV flare of 28 July 2006, whose SED has been 
fitted with r = 50 and with the SSC and EC models (see Fig. 
4). Adapted from [7, 19]. 



3. Matter dominated jets 

Fig. 2 shows the distribution of powers that the jet 
carries in the form of bulk motion of protons (ip, 
assuming one proton per electron), relativistic emit- 
ting electrons (io), magnetic field {Lb) and radia- 
tion (Lr — r^Lj., where L[ is the comoving luminos- 
ity). These powers result from the fitting of a sample 
of blazars for which we have high energy observa- 
tions [7]. Bear in mind a possible bias: we are select- 
ing sources in an active state. All fits assume a sim- 
ple one-zone, leptonic, synchrotron inverse Compton 
model. The different histograms represent FSRQs, 
"classical" BL Lacs and TeV BL Lacs. Note that: 
The jet power, in FSRQs, is large. Estimating 
the luminosity of the accretion disk from the broad 
emission lines, and sometimes from the directly visi- 
ble blue bump emission, we derive that the jet power 
is often greater than the disk luminosity. This is even 
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more true for BL Lac objects, where no lines are seen. 
The power in Poynting flux {Lb) is small. Even 
smaller than L^- This argues against the idea that 
blazar jets are magnetically dominated at all scales 
(e.g. [5, 24]), if the magnetic field of the emitting 
region is representative of the overall magnetic field 
at these scales. A small magnetic field in the emis- 
sion region is expected, since these sources are the 
brightest 7-ray emitters: they emit more by the In- 
verse Compton process than by synchrotron, and this 
limits the value of the magnetic field. 
TeV BL Lacs have the least powerful jets. This 
suggests that it is easier to make them decelerate 
even at the pc scale. 
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Fig. 3. The random Lorentz factor, 7pcaki of the electrons 
emitting at the peaks of the SED vs the comoving energy 
density (magnetic plus radiative; top panel) and the total jet 
power (bottom panel). The star and the triangle correspond 
to the values of PKS 2155-304 during the TeV flare of 28 July 
2006, whose SED has been fitted with F = 50 and with the 
SSC and EC models (see Fig. 4). Adapted from [7]. 



3.1. The blazar sequence 

There are two versions of the so called 
quence' 



blazar se- 

the first is purely phenomenological: the 



observed SEDs show that the two broad peaks (one 
in the IR-soft X-rays, the second in the MeV-GeV 
and sometimes TeV bands) shift to lower frequen- 
cies as the observed bolometric luminosity is in- 
creased. At the same time, the Compton dominance 
(the Lc/Ls ratio) increases [12], see also Maraschi et 
al. in these proceedings). The second version of the 
blazar sequence explains this phenomenological be- 
havior in terms of increased radiative cooling as the 
total luminosity increases. Electrons emitting at the 
peaks of the SEDs have smaller energies as the lu- 
minosity, thus the cooling, increases [14, 16]. Fig. 3 
shows 7peak VS the comoving energy density (mag- 
netic plus radiative; top panel) and vs the total jet 
power (bottom panel) . Different symbols identify FS- 
RQs, "classical" BL Lacs and TeV BL Lacs. The 
general trend is evident: 7poak oc {Ub + U^.)^^ for 
large 7poak, while two branches appear at small val- 
ues of 7pcak [oc (Ub + U'X^ and oc {Ub + U',)-^'\ 
The bottom panel shows another well defined trend: 
7peak oc where = + + Lb- These 

demonstrates (in the context of the used model) a 
tight link between the amount of radiative cooling 
and the power of blazar jets. It also illustrates that 
TeV BL Lacs have the least powerful jets and the 
most energetic electrons. This fact bears an impor- 
tant consequence, as discussed below. 



4. Ultrafast TeV variability 

The increased sensitivity of the new generation of 
Cherenkov telescopes allows to observe flux varia- 
tions down to the minute timescales, for the brightest 
sources. Indeed, PSK 2155-304 [2] and Mkn 501 [3] 
showed variations on ivar =3-5 minutes. Particularly 
challenging is the case of PKS 2155-304, because 
the variations occurred during an overall very active 
state of the source with an observed TeV luminosity 
of ~ 10''^ erg s~^ (see Fig. 4). The usual way to infer 
the size from the variability timescale is i? < ctvar<5 
where 5 is the Doppler factor. Begelman, Fabian & 
Rees [4] pointed out that one consequence of such a 
small ^var is that these timescales are no longer in- 
dicative of the size of the black hole (as instead are 
in the "internal shock scenario" [15, 28]). The other 
consequence is that to avoid to have a too compact 
source, with the accompanying problem of 7-ray ab- 
sorption through the 7-7 — > process [9], one is 
obliged to increase the bulk Lorentz factor of the 
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emitting region (hence 5) to values close to 50 or 
more (close to the ones of Gamma Ray Bursts). In 
turn, such large T makes any external photon source 
strongly boosted in the comoving frame, favoring the 
EC process. At the same time, the viewing angle, to 
have (5 ~ r, must be of the order of 1° or less. There- 
fore Begelman et al. [4] suggest that the EC process 
is the main radiation mechanism. They also argue 
that, since to produce TeV photons one needs highly 
energetic electrons, the jet could be particle starved 
(i.e. one needs fewer electrons to produce the radia- 
tion we see, if they are at high energies). Therefore 
the jet should be magnetically dominated. 
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Fig. 4. The SED of PKS 2155-304. TeV data correspond 
to the flare of 28 July 2006, while X-ray and optical data 
are not strictly simultaneous, but corresponds to 2 days later 
(see [11]). We modelled the SED with an SSC (dashed line) 
and EC (solid line) assuming F = 50, a viewing angle of 1° and 
a size of the emitting region of _R ~ 3 X 10^* cm, correspond- 
ing to fvar ~200 seconds. The dotted line corresponds to the 
distribution of seed photons assumed for the EC model. It has 
been chosen in order to maximize the high energy output of 
the EC emission, without overproducing the far IR observed 
flux. Adapted from [19]. 



Fig. 4 (see [19]) shows the TeV data correspond- 
ing to the flare (observed and de-absorbed, accord- 
ing to a low-medium value of the IR background), 
together to X-ray and optical data taken a few days 
later [11]. The shown models are a pure SSC model, 
with no external photons, and an EC model with 
external photons produced very close to the emit- 



ting region, at frequencies compatible with scatter- 
ings in the Thomson regime, but whose flux does not 
overproduce the observed far IR data of the source. 
In other words, we have maximized and optimized 
the possible contribution of the external component, 
imposing, at the same time, the maximum possible 
magnetic field, to see if the jet of PKS 2155-304 
can be magnetically dominated. For both models, 
r = 50, the viewing angle is 1° and then 6 = 56. The 
size of the emitting region is i? ~ 3 x 10^* cm, yield- 
ing fvar ~ 200 s. Both the SSC and the EC model can 
reproduce the high energy data, but the EC model 
cannot contribute much to the X-ray flux, which can 
instead be reproduced by the SSC model. In Fig. 2 
we show the jet powers of PKS 2155-304 according 
to the two models. Note that Lb is not dominant 
even in the EC model, where the enhanced radiation 
energy density (due to the strong boosting) allows a 
larger magnetic field with respect to the SSC model 
(B = 3 G vs 0.58 G, respectively). Note also that the 
power corresponding to the produced radiation (L^) 
barely corresponds to Lp + for the SSC case, and 
is larger in the EC model. In Fig. 3 we show the value 
7pcak of PKS 2155-304 vs the comoving energy den- 
sities and jet power. The conclusion is that although 
the EC model allows (and requires) a larger mag- 
netic field, the jet remains matter dominated even 
with r = 50. At the same time, the EC model gives 
a worse fit (we cannot account for the X-ray flux) 
with respect to the SSC model, although a variant 
of the EC model, in which the variable TeV emis- 
sion is produced by a very fast "needle" immersed 
in a larger and slower "normal" jet, gives the best 
results [19]. Finally, the SSC model gives values of 
7poak and U' in agreement with the general trend 
observed for other TeV blazars. [7] 

5. Structured and decelerating low 
power jets 

There are several lines of evidence that the jets in low 
power FR I radio-galaxies, hence the jet of BL Lac 
objects, can have a fast inner spine surrounded by 
a slower layer. The evidences come from statistical 
arguments (unifying FR I with BL Lacs [8]), and 
direct radio imaging [21, 22, 25, 31]. Furthermore, at 
the VLBI scale, the jets of low power TeV BL Lacs 
are slow, often subluminal [10, 27]. This is in marked 
contrast with more powerful jets. A slow motion at 
the VLBI scale is also at odd with with the need 
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of bulk Lorentz factors exceeding 20 for the sub-pc 
scale of TeV BL Lacs. At these scales, where most 
of the emission is produced, a large F is required not 
only to avoid the 7-7 process [4, 9], but also 

to account for the large frequency separation of the 
two peaks of the SED [32]. 

To explain the latter properties, Georganopoulos 
& Kazanas [13] have postulated that the jet experi- 
ences a strong deceleration at the sub-pc scale. The 
jet continues to emit during the deceleration phase, 
and the emission at different distances from the black 
hole is beamed differently. The region at the start of 
the deceleration region sees the radiation produced 
farther away along the jet beamed in its direction. 
Viceversa, the regions moving at lower speeds see the 
radiation of the inner jet more beamed. In this way 
there is an enhancement of the radiation energy den- 
sity seen in the comoving frames with a consequent 
enhancement of the inverse Compton emission. 

Ghisellini, Tavecchio & Chiaberge [17], following 
this idea and the above mentioned evidences, pro- 
posed that the jet can be structured, with a cospatial 
fast spine surrounded by a slow layer. Both compo- 
nents emit, and the radiation produced by one is seen 
boosted by the other, enhancing the inverse Comp- 
ton emission of both components. 

This structure influences the jet propagation. 
Since the typical energies of the emitting electrons 
in low power jets are very large (see Fig. 3), the ac- 
companying protons arc not important for the jet 
dynamics. Then if the jet has i) a small power and 
ii) large electron energies, it can decelerate by the 
Compton rocket effect. In fact, the electrons in the 
spine efhciently scatter the radiation produced by the 
layer. In the spine comoving frame, this emission is 
not isotropic (contrary to the synchrotron and SSC 
emission); if the produced radiation is a good frac- 
tion of the electron energy (as it seems, see Fig. 2), 
then the jet decelerates [17]. 

The layer, having a smaller bulk Lorentz factor 
(^a few), emits in a cone much wider than the spine. 
Observers at large viewing angles see the layer, not 
the spine. Since also the layer copiously emits at high 
energies (by scattering spine photons), this idea pre- 
dicts that several radio-galaxies should be detected 
by GLAST [17]. In one case, M87, we might have al- 
ready seen the TcV emission of the layer, as discussed 
below. 
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Fig. 5. SED of the core of M87 together with the H.E.S.S. 
spectra taken in 2004 and 2005 [1]. The upper model is the 
emission from the spine that would be observed at 6° from 
the jet axis, compared with the SED of BL Lac. From [34]. 

5.1. TeV radiogalaxies: the case of M87 

M87 is the only non-blazar extragalactic TeV source 
[1] among the 20 already detected (see www. 
mppmu.mpg.de/'--^rwagner/sources). Even if it is not 
possible to resolve the TeV emission region, the TeV 
flux variability (~ days) suggests a compact source. 
Among the sites proposed as emission regions of the 
TeV flux are the resolved jet (in particular the so- 
called knot HST-1) and the unresolved base of the 
jet, in analogy with blazars. 

HST-1 is located at 60 pc (projected) from the 
core of of M87. It showed spectacular activity in the 
past, and this extreme phenomenology is described 
by [30] assuming that HST-l marks the recoUima- 
tion shock of the jet. As such, HST-1 is thought to 
be a rather efficient particle accelerator and thus a 
possible source of intense TeV radiation, but its lo- 
cation at large distances from the core contrasts with 
the rapid TeV variability, that could instead be eas- 
ily reproduced if the TeV emission is produced in the 
more compact "blazar" -like region. 
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The determination of the TeV emission site is 

rather important: if the emission site will be eventu- 
ally identified as knot HST-1, this will have a broad 
impact on the current view of relativistic jets. 

We [34] have then tried to model the SED of M87 
with a spine-layer jet, assuming that the TeV flux is 
produced by the layer, while the rest of the spectrum 
is produced by the spine, dominating the emission at 
lower energies even at the relatively large M87 view- 
ing angle (~18°, see Fig. 5). The emission region is 
at the same distances than in blazars. "De-coupling" 
the two peaks of the SED (one is made by the spine, 
the other by the layer) solves the most severe diffi- 
culty faced by standard one-zone SSC models when 
applied to the SED of M87. The price to pay is to 
increase (to double) the number of parameters. For 
this reason we not only fitted the M87 SED, but were 
careful to see the model predictions for observers lo- 
cated at smaller angles. As can be seen in Fig. 5, 
the blazar "paired" to M87 (for a viewing angle of 
6°) has a SED closely resembling the one of BL Lac 
itself. This is only a "consistency" check, but it is 
encouraging. 

Acknowledgments 

I gratefully thank Annalisa Celotti, Luigi Foschini 
and Fabrizio Tavecchio. 

References 

[I] F.A. Aharonian, A.G. Akhperjanian, A.R. Bazer- 
Bachi, et al., 2006, Science, 314, 1424 

[2] F.A. Aharonian, A.G. Akhperjanian, A.R. Bazer- 

Bachi, et al., 2007, ApJ, 664, L71 
[3] J. Albert, E. Aliu, H. Andorhub ct al., 2007, ApJ, 

669, 862 

[4] M.C. Begelman, A.G. Fabian & M.J. Rees, 2008, MN- 
RAS, in press (astro-ph/0709.0540) 

[5] R.D. Blandford, 2002, in: "Lighthouses of the Uni- 
verse", ed. M. Gilfanov, R. Sunyaev and E. Ghura- 
zov, Springer-Verlag, p. 381 (astro-ph/0202265) 

[6] A. Bulgarelli, A. Ghen, S. Vercellone, et al., 2007, Atel 
1167 (http:/ /www.astronomerstelegram.org/?read= 
1167) 

[7] A. Gelotti & G. Ghiseilini, 2008, MNRAS, in press 

(astro-ph/0711.4112) 
[8] M. Chiabcrgc, A. Celotti, A. Gapetti & G. Ghiseilini, 

2000, A&A, 358, 104 
[9] L. Dondi & G. Ghiseilini, 1995, MNRAS, 273, 583 
[10] RG. Edwards & B.G. Finer, 2002, ApJ, 579, L70 

[II] L. Foschini, G. Ghiseilini, F. Tavecchio et al., 2007, 
ApJ, 657, L81 



[12] G. Fossati, L. Maraschi, A. Celotti, A. Gomastri & 

G. Ghiseilini, 1998, MNRAS, 299, 433 
[13] M. Georganopoulos & D. Kazanas, 2003, ApJ, 594, 

L27 

[14] G. Ghiseilini, A. Celotti, G. Fossati, L. Maraschi & 

A. Gomastri, 1998, MNRAS, 301, 451 
[15] G. Ghiseilini, 1999, Astronomische Nachrichten, 

320, 232 

[16] G. Ghiseilini, A. Celotti, L. Costamante, 2002, 

A&A, 386, 833 
[17] G. Ghiseilini, F. Tavecchio & M. Chiaberge, 2005, 

A&A, 432, 401 
[18] G. Ghiseilini, L. Foschini, F. Tavecchio & E. Pian, 

2007, MNRAS, 382, L82 
[19] G. Ghiseilini & F. Tavecchio, 2008, subm to MNRAS 

(astro-ph/0801.2569) 
[20] P. Giommi, A.J. Blustin, M. Gapalbi, et al. 2006, 

A&A, 456, 911 
[21] G. Giovannini, G.B. Taylor, E. Arbizzani, 1999, 

ApJ, 522, 101 
[22] M. Giroletti, G. Giovannini, L. Feretti et al., 2004, 

ApJ, 600, 127 
[23] K. Katarzyuski & G. Ghiseilini, 2007, A&A, 463, 

529 

[24] M. Lyutikov & R.D. Blandford, 2003, in "Beaming 

and Jets in Gamma Ray Bursts", Ed. R. Ouyed, p. 

146 (astro-ph/0210671) 
[25] F.N. Owen, P.E. Hardee & T.J. Cornwell, 1989, 

ApJ, 340, 698 
[26] E. Pian, L. Foschini, V. Beckmann et al., 2006, 

A&A, 449, L21 
[27] B.G. Finer & P.G. Edwards, 2004, ApJ, 600, 115 
[28] M. Sikora, M.C. Begelman & M.J. Rees, 1994, ApJ, 

421, 153 

[29] M. Sikora, R. Moderski & G.M. Madejski, 2008, 
ApJ, in press (astro-ph/071 1.3524) 

[30] L. Stawarz, F.A. Aharonian, J. Kataoka, M. Os- 
trowski, A. Siemiginowska & M. Sikora, 2006, MN- 
RAS, 370, 981 

[31] M.R. Swain, A.H. Bridle & S.A. Baum, 1998, ApJ, 
507, L29 

[32] F. Tavecchio, L. Maraschi & G. GhiselUni, 1998, 

ApJ, 509, 608 
[33] F. Tavecchio, L. Maraschi, G. Ghiseilini et al., 2002, 

ApJ, 575, 137 
[34] F. Tavecchio & G. GhiselUni, 2008, MNRAS, in 

press (astro-ph/0801.2569) 



